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Abstract — We investigate the total resistive losses incurred 
in returning a power network of identical generators to a 
synchronous state following a transient stability event or in 
maintaining this state in the presence of persistent stochastic 
disturbances. We formulate this cost as the input-output H 2 
norm of a linear dynamical system with distributed distur- 
bances. We derive an expression for the total resistive losses 
that scales with the size of the network as well as properties 
of the generators and power lines, but is independent of the 
network topology. This topologically invariant scaling of what 
we term the price of synchrony is in contrast to typical power 
system stability notions like rate of convergence or the region 
of attraction for rotor-angle stability. Our result indicates that 
highly connected power networks, whilst desirable for higher 
phase synchrony, do not offer an advantage in terms of the 
total resistive power losses needed to achieve this synchrony. 
Furthermore, if power flow is the mechanism used to achieve 
synchrony in highly-distributed-generation networks, the cost 
increases unboundedly with the number of generators. 

I. Introduction 

Many factors such as increased demand, renewable energy 
mandates [1] and further deregulation of the electric power 
industry [2], [3], [4] are driving changes to the electric power 
grid. The new grid will have to deal with higher levels of 
uncertainty from renewable energy sources and changing 
load patterns as well as increasingly distributed energy 
generation. These changes are likely to make it more prone 
to stability issues. In particular, they have the potential to 
create problems associated with rotor-angle stability, which 
is the ability of the power grid to recover synchrony after a 
disturbance [5], Synchrony refers to the condition when both 
the frequency and phase of all generators within a particular 
power network are aligned. Loss of synchrony can lead to 
load shedding or even black-outs. 

A special case of rotor-angle stability is the so-called 
transient stability problem, which is associated with large 
angle disturbances due to events such as generator, power 
line or other component failures or other abrupt changes that 
can be caused, for example, by intermittent renewable energy 
sources. There is a large body of power system transient 
stability literature, see e.g. [6], [7]. Most of this literature 
focuses on the existence of Lyapunov like energy functions 
[8], [9] and their use in determining a region of attraction 
type criteria for a particular synchronous state or set of states, 
see e.g. [10], [11]. 
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A recent research trend has been to draw connections 
between problems in distributed system control and power 
network stability. This literature is vast, but we note in 
particular a series of works that use a set of coupled Ku- 
ramoto oscillators to represent the power network [12], [13]. 
This non-uniform Kuramoto oscillator modeling framework 
uses a first order approximation of the network reduced 
classical power system model to provide network parameter 
dependent analytical conditions for frequency and phase 
synchronization [13]. Similar first order models have been 
employed to investigate how power flow scheduling and 
adding more power lines (i.e., increasing graph connectivity) 
affects the rate of convergence in a power network [14]. 

In the present paper, we examine the connections between 
distributed control systems and power network stability in a 
different context. We do not study stability, but rather assume 
that the network will return to a synchronized state after 
disturbances. Instead we focus on the cost of keeping the 
network in synchrony, i.e. how much real power is required to 
drive the system to a stable, synchronous operating condition. 
Lack of synchrony leads to circulating currents [15] passing 
between generators whose angles are out of phase. This flow 
of current leads to resistive power losses over the power 
lines due to their non-zero line resistances. This loss is 
generally considered relatively small compared to the total 
real power flow in a typical power network. It is however 
unclear whether these losses will be small in power grids 
of the future, which are expected to have highly distributed 
generation, and consequently many more (though typically 
smaller power capacity) generators than today's grid. 

The problem we analyze is that of a large network of many 
identical generators. We consider several scenarios such 
as the power network encountering disturbance (transient 
stability) events, or being subjected to persistent stochastic 
disturbances where the system is continuously correcting for 
these disturbances. In both of these scenarios, we quantify the 
total power lost due to non-zero line resistances. Our main 
result, in equation ( |13) , shows that this power loss scales 
with the product of the network size and the ratio of line 
resistances to their reactances (which is assumed to be the 
same for all links). The latter quantity is normally assumed 
to be rather small. However, our result shows that even 
though that ratio is small, when the number of generators 
becomes very large, then so will the total resistive losses. 
Furthermore, and perhaps more surprisingly, these losses are 
independent of the network topology, i.e. highly connected 
networks (which would have highly coherent, or close to 
identical, phases) and loosely connected networks (which 



would have a higher degree of phase incoherence) incur the 
same resistive power losses in recovering synchrony. 

Our results indicate that the cost of maintaining synchrony 
using power flows is essentially a function of the number of 
generators in the network and not of its topology. We should 
point out that this conclusion is not inconsistent with other 
results on power system stability and performance measures. 
For example, interaction strength and network topology play 
important roles in determining whether a system can syn- 
chronize [8], [12], [13], [16] and the rate of convergence or 
damping of a power system is directly related to the network 
connectivity [14]. The numerical examples we present in 
this paper confirm that while losses are independent of 
the network connectivity, a highly connected network will 
stabilize more quickly and with less oscillatory behavior. 

One measure that quantifies the lack of synchrony between 
network elements is that of "coherence" [16]. As expected, 
highly connected networks are more coherent than loosely 
connected ones. We refer the reader to [16] for asymptotic 
scalings of lattice-type networks and to [17] for fractal 
networks, and note that these scalings are somewhat different 
than those that quantify network damping. One intuitive 
explanation for the cost of synchrony being independent 
of network topology is as follows. Comparing a highly 
connected network with a loosely connected one, we expect 
the former to have much more phase coherence than the 
latter. Consequently the power flows per link in a highly 
connected network are relatively small, but there are many 
more links than in the loosely connected network. Thus in the 
aggregate, the total power losses of the two networks are the 
same. One should keep in mind however, that a less coherent 
network is more likely to have transient stability problems, 
exit the region of attraction, etc. The issues of stability and 
the cost of synchrony are two different issues. 

The remainder of this paper is organized as follows. 
Section [fl] introduces the classical power system model and 
defines the notation. In Section [ill] we quantify total resistive 
line losses as an input-output H 2 norm of a linear dynamical 
system. This framework is then used to derive an algebraic 
expression for the resistive line losses in terms of the param- 
eters of the admittance matrix, generator damping and the 
size of the network. Various interpretations of the H 2 norm 
are then presented together with the corresponding operating 
scenarios for a power network . The results of the analysis are 
demonstrated for two different network topologies in Section 
IV The paper concludes and suggests directions for future 
work in Section [V] 

II. Problem Formulation 

Consider a network of N buses (nodes) and £ edges. At 
each node i = 1, . . . , N there is a generator Gi, with inertia 
constant Mj, damping ft, voltage magnitude Vi and angle 
8i. Using the classic machine model, which assumes a fixed 
voltage behind a reactance, the dynamics of the i th generator 
is given by [18] 
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Fig. 1 : A depiction of a linear network of n generator nodes 
(buses) Gi for % = 1 . . . , n. 



Here, P e i is the electrical power flowing out of the i th 
generator, and P m .i is the mechanical power input from the 
corresponding turbine. 

Define the admittance over edge £y connecting nodes i 
and j as = gij — j&ij, where g^ and bij are respectively 
the conductance and susceptance of the line defined by the 
edge £ij. Then, the electrical power injection at node i is 
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+ 2>d^l|T4|sin(0i-0fe), (2) 

where k ~ i denotes an edge Eg. and ga is the self 
admittance of the generator. In what follows, we use a 
simplified model in which all generators are assumed to have 
the same self admittance, which we denote g — ga for all 
i — 1, . . . , n. Finally, applying the standard linear power flow 
assumptions used in transient stability analysis (conductances 
are negligible, angle differences are small and voltages are 
constant with unit magnitude) to (|2]i yields 
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Substituting this into ([TJ leads to 
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In order to simplify the notation we define the entries of 



the admittance matrix Y £ 



as 



^Vij, if i = j, 
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if i 7^ j and j ~ i, (i.e., £y 6 £), 
otherwise. 



Prr 



Vi = 1,2. ../j. 



(1) 



Then, we partition this admittance matrix into the real 
(resistive) and imaginary (reactive) parts and define 

Y = Re{F} + jlm{r} =: (L G + gi) + j L B . 

By this construction, Lq and Lb retain the symmetry of Y, 
and they have as a common eigenvector the vector 1 with 
components all equal to 1, i.e. 

L B 1 = L G 1 = 0. 

It is a well-known result that if the graphs underlying the 
system represented by Lb and Lq, which are analogous 



to weighted graph Laplacians, are connected (i.e. any two 
nodes are connected by a path of edges), then the remaining 
eigenvalues of Lb and Lq are all positive. 

An assumption we will invoke later is that all the eigen- 
vectors of Lb and Lq are shared (u is an eigenvector of Lb 
iff it is an eigenvector of Lb)- An important consequence of 
this assumption is that the two Laplacians are simultaneously 
diagonalizable by the same orthogonal transformation U, i.e. 
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where A^ and A^ are diagonal matrices with the respective 
eigenvalues of Lb and Lq as diagonal entries. One setting 
in which this assumption holds is when the ratios of each 
connection's conductance to susceptance (equivalently the 
ratios of resistances to reactances) are all equal, i.e. 

9ik r lk 

— = = a, 

Oik %ik 

which is independent of the link indices (i, k). It then follows 
that 

Lq = aL B , (6) 

which implies that Lq and Lb have the same eigenvectors. 
Finally, we rewrite equation Q in state space form 
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=: Aip + Bw, 



where we have assumed that the constant Mj = 1 for all i 
and that mechanical power input P m .i is a constant that can 
be captured as a part of the input w. 

Remark 1: Based on the classical machine model P mi is 
essentially the real power input at the bus. Thus one can 
consider P m ,i as a system input. Alternatively, the absence 
of P m .i can be justified as follows. For a system with no 
loads the equilibrium point of the system is at P e> i = P m ,i- 
If ([T| defines a linear system (as assumed here) a coordinate 
transform can be applied to obtain a new system with shifted 
angles 9 s corresponding to the transformed system with the 
equilibrium point P m ,i — 0. By abuse of notation we denote 
the shifted coordinates as 9 in the state space representation 
Q. 

Remark 2: The network model |7) assumes a system with 
no loads. Such a model may arise from network reduction 
leading to a system of internal generator buses with loads 
approximated as impedances that have been absorbed into the 
network. The model could easily be extended to include static 
loads through a system of differential algebraic equations 
similar to those derived in [14], [19]. 

III. System Performance 

Several performance metrics can be used to quantify the 
relative stability for the system dTll. A few common metrics 
are the ability to synchronize, the degree of synchronization 
that is achievable, the time to synchronize, or the control 
effort required to obtain the desired system state. In the 
distributed systems setting it is also common to evaluate 
these metrics with respect to various control strategies. For 



example, in evaluating the performance of a distributed 
system with local versus global control strategies. In this 
section we focus on the control actuation required to drive 
a system to the synchronous state. This control effort is 
measured through the real (resistive) power loss over each 
line. These losses are associated with circulating currents 
that arise from the angle differences between generators, i.e. 
disturbances [15]. 

The power flow over an edge £y is 

Pij + Pji = V, {V, - Vj)* Vij + Vj (Vj - VA* y jU 

where * denotes the complex conjugate. The resistive power 
loss over can therefore be defined as 
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Using a small angle approximation and standard trigonomet- 
ric identities this can be approximated as 
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We are interested in the sum total of the resistive losses over 
all links in the network, which is given by 

|2 
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(7) This last quantity can be expressed in vector form as 



Pioss = y*y 

where the vector signal y is an output of the linearized swing 
equations |7]) 
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with CfCi := L G . (11) 



A simple choice for C\ would be to take C\ — L G (possible 
since Lq is positive semi-definite), which is what we assume 
from now on. 

We now calculate the H 2 norm from disturbance w to 
the performance output y of the system |7]i and ( fTTj i. The 
square of the H 2 norm has several standard interpretations 
including (a) The variance of the output y when the input 
w is a unit variance white stochastic process, (b) The total 
time integral of the variance of y when the initial condition 
is a random variable with correlation matrix BB* , and (c) 
The total sum of time integrals of output response powers 
given an impulse as a disturbance input at each generator. 

We first perform the H 2 norm calculation and derive a 
formula in terms of the system's parameters, and then inves- 
tigate the implications of the three different interpretations 
for this particular system of swing equations. 

H 2 Norm Calculation 

For ease of reference we rewrite the system equations |7]) 
and ( fTT| here 
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We will denote the input-output mapping of this system 
by H. This system has all eigenvalues strictly in the left 
half of the complex plane with the exception of a single 
zero eigenvalue of Lb- It is, however easy to show that 
this unstable mode is unobservable from the output y due 
to the fact that Lq = C\C\ shares this eigenvalue and 
its corresponding eigenvector (see Appendix for the full 
argument). It then follows that the input-output transfer 
function from w to y is indeed stable and has finite H 2 
norm. 

The easiest method to calculate the H norm of the <[12J» 
is using a spectral decomposition of L B - Consider the state 
transformation 
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where U is the matrix in |5]l diagonalizing Lb and Lq. The 
system in the new state variables becomes 
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Since multiplying by orthogonal matrices does not change 
the H 2 norm, we can multiply the input by U and and the 
output by U* (i.e. define w = U*w and y = U*y) to obtain 
an equivalent system (that has the same H 2 norm) 
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We will denote the input-output mapping of this system 
by H. Since A# and Aq are diagonal, this represents N 
decoupled systems 
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where n = 1, . . . , N are the indices of eigenvalues A,^ and 
A„, which correspond to Lb and Lq respectively. Denote 
the input-output mapping of each decoupled subsystem by 
H n ■ We can then write 



H = diag H l 
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The square of the H 2 norm of the system ( [T2] i is thus the 
sum of the squares of the H 2 norms of all the decoupled 
subsystems 
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The H 2 norm of each of the individual subsystems can now 
be easily calculated as follows. For n = 1, the corresponding 
eigenvalues are Af = Af = 0, and we have a completely 



unobservable system, thus H-H1II2 = 0. For n 7^ 1, we solve 
the Lyapunov equation for the observability Grammians X n 
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where for simplicity of notation we have dropped the sub- 
script n from the components of of the Grammian X n . This 
matrix equation corresponds three equations, of which only 
the following two are relevant 
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Finally, since the B matrix of each subsystem is [0 1] , the 
H 2 norm (squared) of each subsystem is just X22, and we 
therefore conclude that 
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The total H 2 norm (squared) of the overall system ( |T2] > is 
thus 
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In summary, we conclude that the total resistive losses in 
this network are 
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where /3 is a generator's self damping, ^ is the ratio of a 
line's resistance to its reactance (assumed equal for all lines), 
and N is the number of generators in the network. Note that 
this expression is independent of the network topology but 
depends on the number of generators in the network. 

H 2 Norm Interpretations for Swing Dynamics 

In the previous subsection we calculated the H 2 norm of 
the linearized swing dynamics |7]), together with the output 
equation (jTTJ based on a disturbance input (forcing) w. In 
our formulation the square of the Euclidean norm y*y of the 
output vector was defined to be equal to the total real power 
dissipated due to resistive losses in the network connections 
as the system synchronized. This synchronized state (the 
equilibrium) is maintained by circulating currents leading to 
power flowing back and forth between generator nodes. 

It is well known that the H 2 norm has at least three 
different interpretations. We recall these interpretations here 
in order to give three different physical scenarios in which 




(a) 5 Node Line Graph (b) Fully Connected 5 Node Graph 





(c) 20 Node Line Graph 
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Fig. 2: (a),(b) Simulation of a 5 bus system of identical oscillators with a random initial velocity disturbance, and zero 
initial phase disturbance. (c),(d) Simulation of a system with 20 identical oscillators under the same conditions as (a),(b). 
Figures (a),(c) are for linearly connected networks (with line graphs as in Figure [TJ, while Figures (b),(d) are for fully 
connected networks (complete graphs). Note how linear networks are much less coherent than complete networks (Figures 
(a),(c) versus (b),(d) respectively). Total resistive losses are however the same for the 5 node networks (a) and (b), and for 
the 20 node networks, (c) and (d), respectively, reflecting that resistive line losses are a function of the network's size rather 
than its topology. 



equation ( fT3| quantifies the resistive losses due to the syn- 
chrony requirement. Denote by H the following Linear Time 
Invariant (LTI) system 

ip(t) = Atp(t) + Bw(t) 

y (t) = cm, 

The H 2 norm ||if||2 of the system has the following three 
different interpretations. 

(a) Response to a white stochastic input. When the input 
w is a white second order process with unit covariance 
(i.e. E {w(t)w*(t)} = /), the H 2 norm (squared) of the 
system is the steady-state total variance of all the output 
components, i.e. 

\\H\\ 2 = lim E{y*(t)y(t)}. 

t— f oo 

For the swing dynamics ([12} the disturbance vector 
can be thought of as persistent stochastic disturbances 
(forcing) at each generator. These disturbances, which 



are uncorrelated across generators, can be due to un- 
certainties in local generator conditions, sudden changes 
in load, or fault events. The variance of the output is 
exactly the expectation of the total power loss due to 
line resistances, 
(b) Response to a random initial condition. When the input 
is turned off, but the initial condition is a random variable 
x a with correlation E {x x*} = BB* , then the H 2 norm 
is the time integral 

/>oo 

\\H\\l = / E{y*(t)y(t)} dt 
Jo 

of the resulting response y. 

The interpretation for ( p~2] > is as follows. Since BB* = 

^ , the corresponding initial condition corresponds 

to each generator having a random initial velocity per- 
turbation that is uncorrelated across generators and zero 
initial phase perturbation. Then \\H ||| quantifies the total 



(over all time) expected resistive power losses due to the 
system returning to a synchronized state, 
(c) Sum of responses to impulses at all inputs. Let 6j refer to 
the vector with all components zero expect for 1 in the 
i th component. Consider N experiments where in each, 
the system is fed an impulse at the i th input channel, 
i.e. Wj(t) = Cj<5(t). Denote the corresponding output 
by yi. The H 2 norm (squared) is then the sum total of 
the L 2 norms of these outputs, i.e. 

\\H\\l = E / #(*)!«(*) dt. 

A stochastic version of this scenario corresponds to 
a system where the inputs w; can occur with equal 
probability. Under this assumption the \\H\\?, becomes 
the expected total power loss given these inputs. 
The interpretation for ( fj"2| ) is when each of the generators 
are subject to impulse force disturbances (since enters 
into the momentum equation of each generator), and 
|| H || I is thus the total power loss over all time under 
such a scenario. 

IV. Numerical Examples 

Consider two networks of identical generators, one whose 
underlying graph is a line, as in Figure [T] and one with 
a fully connected graph. In this section, we compare the 
behavior of systems with these two network topologies as 
the network size varies. All simulations use the following 
parameter values [20]: M — ^K, (3 — with a frequency 
/ = 60 Hz. The admittance between connected generators 
is set to yij = 0.2 + j 1.5 and we assume g it =: g = for 
all i = 1 n. 



Figures 2a and 2b respectively show the state trajecto- 
ries of a 5 node system with a line graph and one with 
a fully connected graph given identical initial conditions. 
The initial conditions for each Wi(0) were drawn from a 
uniform distribution [—1,1], and 0(0) was set to zero. This 
corresponds to the H 2 norm interpretation (b) described 
in the previous section. For both the line graph and the 
fully connected network the total losses are the same, with 
pioss _ i 2885, as predicted by the main result in ( fH} . 
However, the transient behavior shows that the system with 
the fully connected graph has more "coherent" phases (i.e. 
they stay closer together). This is most clearly visible in 
the top panels of the plots that depict the state 0. Here, the 
small oscillations for the line graph system continue through 



the 10 second interval shown in Figure 2a but have almost 
completely died out for the fully connected of Figure 2b 

In order to evaluate the effect of increasing the network 
size we ran a similar test with two 20 bus systems with the 
same topological structure and initial conditions (u>i(0) G 
[—1,1] drawn from a uniform distribution for each i and 
0(0) = 0). In this case, the losses increased to p loss — 
5.9256 but remained equal in both networks. The faster 
convergence to the synchronized state in fully connected 
system (compared to the line graph) is much more evident 



The simulation results indicate that a fully connected 
graph is much more coherent than a line graph. The resem- 
blance of this result to those about coherence in vehicular 
formations (platoons) and similar consensus-like network 
algorithms [16] is striking. With coherence as a performance 
measure, the line graph is the worst such topology (amongst 
connected graphs), while the fully connected topology is the 
best. However this additional coherence, while desirable for 
other reasons, will not result in lower resistive losses. 

V. Conclusions 

We have considered a power network model with dis- 
tributed disturbances, and quantified the total resistive power 
losses incurred due to the current flows needed to maintain 
phase synchrony. We have shown that these losses are in- 
dependent of the network's topology, but scale unboundedly 
with the number of generators. There are interesting impli- 
cations for the design of future highly-distributed-generation 
networks, which have potentially orders of magnitude more 
generator nodes that today's networks. While a highly phase 
coherent (and therefore highly connected) network is desir- 
able for many reasons, the cost of maintaining this coherence 
depends only on the number of generators and not the 
network's connectivity. Since this cost grows unboundedly 
with the number of generator nodes, then the current scheme 
of using power flows as the synchronization mechanism may 
not be scalable to future networks. This is perhaps a further 
argument for the use of other control mechanisms, such as 
communication links, for phase synchronization. 

Appendix 

Proof that \\2\ is input-output stable 

It is well known that for any pair (C, A), the observability 
of (C, A) and (C*C,A) are equivalent. The only unstable 
mode of the A-matrix in ( fL2j i is at zero with corresponding 



eigenvector 



This mode is unobservable from the output 
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in the larger network, as shown in figures 2c and 2d 



y since by the PBH test 
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Note that this is a consequence of Lq and Lb having the 
common eigenvector 1 with eigenvalue 0. The unobservabil- 
ity of the only unstable mode implies then that the system 
( fT2| > is input-output stable. 
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